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Abstract. In this paper, we deal with nonlinear ill-posed problems involving 
monotone operators and consider Lavrentiev’s regularization method. This ap¬ 
proach, in contrast to Tikhonov’s regularization method, does not make use of 
the adjoint of the derivative. There are plenty of qualitative and quantitative 
convergence results in the literature, both in Hilbert and Banach spaces. Our 
aim here is mainly to contribute to convergence rates results in Hilbert spaces 
based on some types of error estimates derived under various source condi¬ 
tions and to interpret them in some settings. In particular, we propose and 
investigate new variational source conditions adapted to these Lavrentiev-type 
techniques. Another focus of this paper is to exploit the concept of approximate 
source conditions. 


1. Introduction. For an infinite dimensional and separable real Hilbert space X 
with norm || • || and inner product (•,•), we consider the ill-posed operator equation 

F{x) = y (1) 

with the nonlinear forward operator F : R(F) C A —>■ A and assume that only 
noisy data G X are available such that the deterministic noise model 

Il2/-/ll<^ (2) 

with noise level <5 > 0 applies. In this context, let G R(F') denote an exact 
solution to equation (1) and x G X a reference element (initial guess). Throughout 
this work we suppose that E is a monotone operator, i.e. we have 

(F(x) — F(x),x — x) > 0 for all x,xGT>(F). (3) 

Then the operator equation (1) is well-posed if 

(F(x) — F{x),x — x) > ||a; — i|| 0(||a: — i||) for all x,x G 'D{F) 
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holds for some index function 0. A function 9 : (0, oo) —>■ (0, oo) is called index func¬ 
tion if it is continuous, strictly increasing and moreover satisfies the limit condition 
limt^+o 0{t) = 0. 

If well-posedness of (1) fails, a regularization approach is required in order to 
find stable approximate solutions to the ill-posed equation (1). We are going to 
construct regularized solutions to by solving the equation 

F{xi)+a{xi-x)=y^. (4) 

This approach, which was originally suggested in [27], is often called Lavrentiev’s 
regularization method. Alternatively, it is called perturbation method (cf. [28]) in 
order to refer to the manner of obtaining regularized solutions as singular pertur¬ 
bations. It can be shown that equation (4) has a unique solution xf^, under various 
types of assumptions. We refer below to two variants of such assumptions, namely 
(i) in the context of a globally monotone operator F defined on the whole space X 
and (ii) when the domain of F might be smaller than X and F is locally monotone 
there. 

(i) If F : A —?► A is globally monotone, i.e. (3) is valid with 27(F) = A, and F is 
a hemicontinuous operator then, by the Browder-Minty theorem, (4) has a unique 
solution x^ for all x S A and € X. 

(ii) If there is a ball Br{x^) C 27(F) around a solution x^ to equation (1) with 

radius r = ~ ^11 + a ^ C A —>■ A is a Frechet differentiable operator 

monotone and hence hemicontinuous operator in the ball, then (4) has a unique 
solution in for all a; G A and y^ € X (cf. [37, Theorem 1.1]). 

In both cases, xf^ depends continuously on y^. 

Note that we will work under assumptions (i), for simplicity of exposition. A 
comprehensive study of Lavrentiev’s regularization method for nonlinear equations 
in Hilbert spaces with monotone operators, even in a more general setting, can 
be found in the book [2], for modifications of the method see also [8, 9, 14, 25]. 
Nevertheless, for completeness of exposition and since some of the estimates will be 
needed later on, in Section 2 we provide a short summary of the arguments leading 
to convergence of from (4) to x^ (see also [1, 6]). 

As in the case of Tikhonov’s regularization method (cf., e.g., [35, Section 3.2] and 
[20]), also for Lavrentiev’s regularization method a certain additional smoothness 
of x^ with respect to the forward operator is required in order to derive convergence 
rates. At the origin of such studies, source conditions (range conditions) of Holder- 
type (cf., e.g., [25, 28, 37, 38]) were under consideration, however attaining here the 
specific form 

x^ — X €'1Z{A^), 0<p<l. (5) 

Extensions of type 

x"^ — X & 'R.{ip{A)), (6) 

with index functions </?, can also be found in the literature (cf., e.g., [29] and [5, 36]). 
In this context, A G £(A, A) denotes a bounded linear operator related to some 
derivative F'{x^) of the forward operator at the solution point. If these conditions 
exceed the present solution smoothness, some weaker forms of source conditions 
might be of interest. 

Recent developments regarding convergence rates theory involve variational in¬ 
equalities as source conditions for Tikhonov regularization, as introduced by [22]. 
One focus of this work is to investigate in case of Lavrentiev regularization not only 
the already known conditions in a Tikhonov-type context, but also new conditions 
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adapted to the Lavrentiev’s techniques - see Section 3 in this respect. Another focus 
is to exploit the concept of approximate source conditions in the sense of [17, 19] in 
Section 4. Finally, we provide some illustrative examples in Section 5. 


2. Convergence. Similarly to the Tikhonov regularization for nonlinear operator 
equations in a Hilbert space setting (cf. [12, Sect. 10.2], also for Lavrentiev regu¬ 
larized solutions convergence is global, in the sense that no sufficient closeness 
between x and needs to be required. Also a general condition which is well- 
known from Tikhonov regularization occurs again here, namely we assume that F 
is weak-to-norm sequentially closed, which means that 

Xn ^ X and F[xn) —>■ => i G F>{F) and F(x) = zq. (7) 

Monotonicity and hemicontinuity of F are together sufficient for (7) in the case 
F>{F) = X (cf., e.g.. Lemma 1.4.5 and Theorem 1.4.6 in [2]). 

We assume that the set L of solutions to equation (1) for given y G X is not 
empty. We derive below formulas for handling the error of regularized solutions x^ 
with respect to exact solutions x^ G L in order to show convergence for appropriate 
choices of the regularization parameter a, where in particular ai-minimum-norm 
solutions are targeted. As usual, we say that a solution G L is the s-minimum- 
norm solution if [[xl^j — a;|| = min — ijj. 

x+ei 

Under condition (7), the set L is weakly closed and hence such a minimum- 
norm solution exists (see, e.g., [35, Prop. 3.14]). If in addition L is convex and 
closed, then is uniquely determined. This is actually the case throughout this 
paper, as the context (i) that we work with, namely F>{F) = X, monotonicity and 
hemicontinuity of the operator, implies maximal monotonicity of F and is sufficient 
for the convexity and closedness of L (cf. Th. 1.4.6 and Corollary 1.4.10 in [2]). 

Testing (4) with x^ — x^ and moreover with F{x^) — F{x'^) leads to 

{F{xi) - F(a;1'), xi - a;^') + {y - y\ xi - 

Fa\\xl^-x'^W'^ + a{x^ - x,xl^-x"^) = Q ( 8 ) 


and 


\\F{xi) - F(xt)||2 + {y- y\ F{xi) - F{x^)) 

+a{F[x^^) - F{x^), - a;’'’) -I- a{x^ - x, F(a;^) - ^’(a;''')) = 0 , 


respectively. By using the monotonicity condition (3) of F in combination with the 
Cauchy-Schwarz inequality this implies the following three basic estimates, which 
will be required below. 


Wxi-x'ff < (x^ - x,x1' -x^) -b -||x^ -x+ll, 

a 

-S|| + 

a 

|]U(xl)-F(xt)||<a||xt-x||+<5. 


(9) 

( 10 ) 

( 11 ) 


Taking into account these estimates one can conclude convergence with a priori or 
a posteriori parameter choice by standard arguments of regularization theory. 
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A priori parameter choice. We choose a» = a* (<5) a priori such that it satisfies 
the limit conditions 

a* ((5) —?► 0 and ——>■ 0 as d —?► 0. (12) 

Q!*(d) 

Inserting this into (10) we obtain that ||a;^ — a;^|| is uniformly bounded as d —>■ 

0. Thus, for any sequence —>■ 0 there exists a weakly convergent subsequence 

weak limit x. On the other hand, (11) and (12) yield V 

and hence, from the weak sequential closedness (7), we conclude that x solves (1). 
Thus, we have weak subsequential convergence to a solution. Convergence is even 
strong since by (9), (10), which remain valid for any solution x of (1), we have 


limsup \\x^Z\sr.^)~ 

k—^oo ^ 


< limsup 

k—^oo 


{x — : 




) + 


C^*(<^nfc) 



= 0 , 


where we have used weak convergence, (12) and uniform boundedness of ^ — 

i||. Along the lines of the proof of [2, Theorem 2.1.2], we show that every such limit 
X is the x-minimum-norm solution. Let x^ be an arbitrary element of L. From (9) 
we have 


(x^ — 


X, X' — X 




and hence, under (12) with x 


Sr 

cx 




—^-xGLasfc—>- 00 , the following holds: 


(x^ — x,x^ — x) >0 for all x^ G L. 

Due to the convexity of L, one has Xt := tx + (1 — t)xl^^^ G L, for all 0 < t < 1. By 
setting x^ := Xt 

(xt - X, (1 - t)(x|„j„ - x)) > 0 for all 0 < t < 1, 


one has 

(xt - X, xl^i^ - x) > 0 
which yields, for the limit t —> 1, 

(x-x,x|„„-x) > 0 and (x - x, - x) > ||x - xf. 

This, however, provides us with the inequality ||x — x|| < ||xj,^j„ — x|| which shows 
that X is the x-minimum-norm solution. By a subsequence-subsequence argument, 
the whole sequence converges to The discussion above proves the following 

proposition. 

Proposition 1. For an a priori parameter choice a, = a* (5) satisfying (12) and a 
sequence i5„ —>■ 0, the sequence of associated regularized solutions x°ff ^ converges 
strongly to the x-minimum-norm solution x|„j„ of (1). 
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A posteriori parameter choice: a variant of discrepancy principle. By in¬ 
spection of the detailed considerations in the previous paragraph concerning the 
convergence for a priori choices of the regularization, we immediately find the fol¬ 
lowing extension to a posteriori parameter choices. 


Proposition 2. For an a posteriori parameter choiee a* = well-defined 

for sufficiently small <5 > 0 and satisfying the conditions 

c 

— 77 —7T- —>■ 0 as S 0 (13) 

a4S,y^) 


as well as 

ll-^(a^l( 5 )) - y|l ^ 0 as 5^0, (14) 

the sequence of regularized solutions ^ associated with a sequence > 0 

converges strongly to the x-minimum-norm solution of (1). 


Now we are going to apply Proposition 2 to a variant of the discrepancy principle 
(compare to the residual principle in Chapter 3 of [2]) and assume in this context 
that 

||F(a:) — y\\ > 0, i.e. x does not solve (1). (15) 

Given a sufficiently small noise level <5 > 0 and data y^ G X satisfying (2), we fix 
0<K<1, T>1 and Oq > 0 such that 

linO-/ll>r5^ (16) 

Under (15), the existence of such an oq for sufficiently small <5 > 0 is due to the 
limit condition 

lim \\F{xi) - y^W = ||P(x) - y^|| 

ct—^oo 

and ||U(a:) — ?/‘^|| > ||i^(a:) — y|| — 5 > for S < min(l, j ). Moreover, 

we know that there is some a > 0 such that 

~ y^W — 0 < a < a. 

This is a consequence of (4), (10) and the estimate 

- /II = a\\xi - x|| < a (11/ - ill -f ||xl - i||) < 2a\\x^ - i|| -f (5, 


which yields a = 

In the following we choose a* = afi5,y^) along a fixed geometrically decaying 
sequence ak = q^ao for some q G (0,1) according to 

a* = maxjafe > 0 : ||P(/j - /|| < rd'^}. (17) 

This is a variant of the sequential discrepancy principle, see also [4]. Well-definedness 
and strict positivity of a* follow from the existence of positive values ag and a. Es¬ 
timate (11) together with (17) yields, for 5 < 1, 

< ||F(a:k)-/|| < ^ ||a;t _ ^|| + j ^ ||^t _ ^|| + 

9 q q 


and thus, the following lower bound for a* 

q{T - 1) 


a* > 


II.tI — 


S^. 
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Consequently, for sufficiently small i5 > 0 one has 


JC < -^11 

a* “ q{T - 1) 


0 as 6 —y 0, 


implying (13). To see (14) from (17), we use the triangle inequality as 

l|i"(^l)- 2 /ll< ll^’(a^ij-/ll+<5<(T + l)<5'‘^0 as ^^0. 


Hence, Proposition 2 applies for this specific variety of a discrepancy principle and 
arbitrary choices of the exponent 0 < k < 1. Note that one can show only weak 
convergence ^ ^ ^Inin as n —?► c» if the exponent is k = 1. 

As an alternative a posteriori choice of the regularization parameter for Lavren¬ 
tiev’s method we can also exploit the Lepskii or balancing principle (cf., e.g., [30]), 
which even yields convergence rates (see Theorem 5). 


3. Convergence rates under variational source conditions. As mentioned 
in the introduction, we investigate below error estimations under several types of 
source conditions on solutions of (1). Actually, the solution referred to in all these 
source conditions has to be the ai-minimum-norm solution in the current setting of 
globally monotone operators F defined on Hilbert spaces. However, we prefer the 
notation in order to anticipate the more general situations when the solution set 
L of (1) is not necessarilty convex and closed. 


3.1. Variational source condition from Tikhonov regularization. First we 
consider the variational source condition 


(x'^ — x,a;^ — x) < f3\\x^ — x\\^ + tl){\\F{x) — F{x 


t'l 


for all X G M, (18) 


which up to some factor 2 on the left-hand side coincides (for ^(t) = t) with the 
variational source condition originally developed in [22] for obtaining convergence 
rates in Tikhonov regularization. Here we use a constant j3 G [0,1), some index 
function ijj and a solution of (1). The set A1 C V in (18) must contain all 
regularized solutions x^ of interest for <5 > 0 sufficiently small and a > 0 chosen 
appropriately. For instance, the results in Section 2 guarantee that, under condition 
(7) and with an appropriate choice of a in dependence of <5 and a fixed but possibly 
small radius p > 0, the set M = Bp{x^) contains all x^ for <5 > 0 sufficiently small. 
By homogeneity argument in the case of a linear operator F = A G C{X,X), we 
conclude that (18) is realistic only for index functions decaying to zero not faster 
than ■!/)(t) ^ t as t ^ -bO. Precisely, taking x := x'^' -b e(x — x'^’) with e > 0 in (18) 
implies 


Ix’^ — x"^ 


11^ < — a^ll^ + 4’{e\\A{x^ — x) 


i.e., unless the trivial case x''' — x G A/’(A) holds, we can divide hy t ■= e [j A(x'l’ — x) 


to conclude that 


for r — 

e - ||A(a;T-S)|| 


lim 

t —^+0 


i’it) 


> c> 0 


(19) 


> 0 . 


As a matter of fact, if 0 < /3 < ^ and A4 contains an x-minimum-norm solution, 
then all solutions x^ satisfying (18) have to be x- minimum-norm solutions. Namely, 
using (18) with x = implies 


1 


“2 - 


- t''' 

min 

i||x'^-x|P- 


1 


2 11^ ^min I 


= (x'^ —i 


t t 

r' — nr'' 


^)< 
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i.e., with 0 < /3 < we get ||a;t - xW^ < - xf. 

Theorem 3. Let F : X ^ X be a hemieontinuous and monotone operator and A4 
he a subset of the Hilbert space X such that (18) is satisfied. Moreover, let a he 
chosen such that all regularized solutions x^ from (4) belong to A4 for sufficiently 
small S > 0. Then one has the estimate 

Iki+ Y^V'(<^ + ca), (20) 

with c = ||a;^ — a;|| for S > 0 as above. For the a priori choice 

a(5)~$-i(52) (21) 

with the index function $( 0 ;) := a^'ip{a) which satisfies (12), this yields the conver¬ 
gence rate 



The best possible rate under condition (18) occurs for tf{t) ~ t as 

IK^)-^^II=0(<5^) */ a(^)^<5i (23) 

Proof. Let 7 := {F{x^) — F{x^),x^ — a;l). Note that x^ G M for sufficiently small 
(5 > 0, whenever the a priori parameter choice a = a{S) satisfies condition ( 12 ) 
guaranteeing (due to Section 2) convergence of x^ to x^ and if A4 = Bp{x^) for 
some p > 0 . 

From ( 8 ), (18) and (2) it follows that 

7 + a|ki= {y^ - y,xi, - x"^) + a{x'^ - x,x'^ - xi) (24) 

< <5||xi - xt|| + a{fi ||xi - xtf + f{\\F{xi) - F{x^)\\)). 

The monotonicity (3) of F implies 7 > 0 and thus we can neglect 7 and estimate 
further as 

(1 - 13) ||xi - xt f < -||xi - xt|| + fiimxi) - F(xt)ll). (25) 

a 

By monotonicity of 'i/? and (11) this yields 

Iki-x+f < -^-\\xi-x^\\ + -^f{5 + ca). 

1 — p a 1 — P 

Since for any d,Ci,C 2 > 0 

d^ < 2cid + C 2 => d^ < Ac\ + 2 c 2 (26) 

holds, one obtains (20). Since $ is an index function, the same is valid also for 
as well as for the compositions $“^((5^) and y/ Consequently, 
the a priori regularization parameter choice ( 21 ) satisfies S = 0(a(S)) and moreover 
the limit conditions (12), which are relevant for convergence (cf. Section 2). Then 
the convergence rate (22) can be derived from (20) for the choice (21). Due to (19), 
we have the inequality < ci$(a), hence $“^((5^) < and , 

for some positive constants c \, C 2 and C 3 . This shows that 0(p-^^) is the best possible 
rate attained for iff) ^ t. □ 
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Remark 4. The classical ‘adjoint’ source condition — x G TZ{A*) (which implies 
(18) with ^ = 0 and ipit) = t in the linear operator case) yields the convergence 
rate Thus, the rate mentioned with no proof by Theorem 6 in 

[28] seems to be a typo. However, the ‘no adjoint’ source condition x"'^ — x G T^{A) 
occurring in (5) for p = 1 implies the better rate (cf. [37]). This is due to 

the higher rate 0{a) for the approximation error (the exact data case), as compared 
to derived from x^ — x G TZ{A*) (see, e.g.. Theorem 5 in [28]), while the 

propagated data error is the same in both situations: 0(5/a). 


3.2. Variational source condition specific to Lavrentiev regularization. In 

this section we propose a variational source condition which seems to naturally fit 
the setting of Lavrentiev’s method. 

Recently, results in [3] regarding Tikhonov regularization in the linear case point 
out connections between the standard source condition G TZ((A*A)''^^) for iz G 
( 0 , 1 ] and several types of variational source conditions along with the one 

2(a;^,a;) < ,5il|a;l|^ +/ 32 |[ 7 la;l|^ for all x G X, (27) 

with constants /3i G [0,1) and /?2 > 0, which is very similar to (18) for power-type 
functions il). Searching for an adaption to the Lavrentiev setting, it makes sense 
however to consider (28) instead of (27), 

(x^,x) </3il|xl|^-I-/ 32 (^a:,a:)'^ for all x G X, (28) 

with fi > 0 and constants Pi, 132 as above. Actually, we will deal with the more 
general version with respect to the nonlinear case 

{x^^ — x,x'^ — x) < P\\x —x'^W'^+ ip{{F{x) —F{x'^),x —x"^)) for all x G M, (29) 

with P G [0,1) and (p an index function. Again, no condition additional to mono¬ 
tonicity needs to be imposed on F. As in Section 3.1, there is a natural upper 
bound on the decay rate of p as t tends to zero. To this end, consider again the 
linear case F = A G C(X,X) and x •.= x^ + £{x — a:'!') with £ > 0 in (29). This 
implies 

< /3£^ Jja:’'’ — Sjp -|- ip{£^{A{x'^ — x),x'^ — x)), 
which by division to i/t = (A(a;f — x), x^ — x) and by letting e —>• 0 enforces 

lim > c > 0 (30) 

t^o - 

with the constant c = , > 0 , unless the trivial case 

yJ{A(x^-x),x\-x) 

{A{x^ — x),x^ — x) = 0 holds, so that the best rate obtainable here is the one 
associated with p(t) = \/t. 

To quantify convergence rates, we introduce the index function 41 as follows. 
Let / be an index function such that its antiderivative /(s) := f{t)dt satisfies 
the condition /(p(s)) < s for s > 0 and let G{a) > /“/“^(T)dT. Then we set 

We will first of all obtain rates by choosing 

Ctapri = cx{S) := 0“^((5^) (31) 


with 


0 (a) := 0^41(0) 


( 32 ) 
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which also satisfies condition (12) for convergence without source condition. 

The same rates are achieved with an a posteriori parameter choice aLep according 
to the Lepskii principle, which does not require the knowledge of T. In this context, 
the regularization parameter aLep is determined as follows: For prescribed r > 1 
and 0 < g < 1 we consider the function 


V3-2f3 \ d 
I- 13 ) a 

and the increasing geometric sequence 


S](q;, S) := 


Ag := {aj : Uj = j G N} 

with ao > 0 sufficiently small, such that \\x^ao ~ ^^^11 — S(ao,(5). Then we set 

otLep = 0t{5,y^) := 

maxja G Ag : ||x^/ < 2 E(q:',( 5) for all a' G A, fl [aojCi]}- 


Theorem 5. Let F : X ^ X be a hemicontinuous and monotone operator and Xi 
be a subset of the Hilbert space X such that (29) is satisfied. Moreover, let a be 
chosen such that all regularized solutions from (4) belong to Xi for sufficiently 
small (5 > 0. Then one has the estimate 

11-^-'It s (3^) 

for 5 > 0 as above. For the a priori choice (31), this yields the convergence rate 

IkL... - II = o ( = o (\/^(e^W), (35) 

with 0 as in (32). The best possible rate occurs for ip{t) ^ and 4'(t) ^ t as 

lKp..-^^ll=0(<5^) */ a{S)^Si. (36) 

The convergence rate (35) is also obtained with the Lepskii principle (33) 

IkLp-^^ll = • 

Proof. Note that G A4, for all 5 > 0 sufficiently small, due to the convergence 
results in the a priori parameter choice case. Let again 7 := {F{x^^) — F(a;i),x^ — 
xi) > 0. By (24) and (29) written for x := x^ one has 

7 + a||xf, -xif < (5||xf, -xi|| + a(/3||xf, -x^f + (^( 7 )). 


Now one can employ Young’s inequality 

pa nb 

ab< f{t)dt+ / f~^{t)dt 
Jo Jo 

for an index function / which will be specified later, for a = tpfj) and b = a. Using 
the settings 

/(s) = f f{t)dt, G{a)> [ f~^{t)dt, 

Jo Jo 

one needs to find the function / such that f{(f{'y)) < 7. Then one has 
a(l -/3)l|a:a < d\\xi-x^\\+G{a), 
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so with 'I' (a) = one gets 

We use again the implication (26) to conclude (34). The a priori parameter choice 
(31) yields the convergence rate (35). For showing the same rate for the a posteriori 
parameter choice according to the Lepskh principle, we apply Proposition 1 from 
[30] (see also [23, § 4.2.2]). Taking into account that dt is an increasing function, 
we have from (34) that 


II 5 t||2^ +2(1 f])ag^pj,^"i'{aapri) _ + 2(1 —/3)i5^ 

- (l-/3)2a2 (l-/3)2a2 

for all ao < a < aapri- This ensures the estimate — a;^|| 

O (cf. [30, Proposition 1]). 


= (E(a,^))2, 

+ 3 ^(^Ctaprif 


□ 


Remark 6 . The construction of f,G and T in Theorem 5 is straightforward if 
ip is smooth enough and we can use f{t) := {if~^y{t) and G{a) := /p“/“^(r)dr. 
Obviously, 'k in Theorem 5 plays a similar role as '0 did in Theorem 3. 


3.3. The special case of Holder and logarithmic source conditions. 

Holder-type variational source conditions. Taking into account exponents /i G (0, i], 
we consider (29) with 

fit) = , f~^it) = /(s) = sf, 

G(q;) = (1 —/r)/xT^aT+‘, 'i'(a) = {1 — p.)^ Q(t) ^ t i-j* , 

or correspondingly (18) with 

0 (t) ^ . 

According to Theorems 3 and 5 this yields in both cases the Holder rate 
\\xi^s)-x^\\=0{Sx^) if a(<5)^<5^. 


Logarithmic-type variational source conditions. We consider (29) with 

= /W=e-4e-=^-'(t), = 


G{a) = 


1 


— In a 


> 


dt, 0 (t) ~ 


— Int 


— Int 


for a,t G ( 0 , 1 ) or correspondingly (18) with 

ipit) 


— Inf 

According to Theorems 3 and 5, this yields again in both cases the same logarithmic 
rate O ■ 

Since the a priori choice a(5) ~ 0 “^(^ 2 ) can not be determined explicitly in this 
logarithmic case, a more convenient choice is a((5) ~ \f5 which implies 


For other results concerning logarithmic rates based on variational source conditions 
we also refer, for example, to [10, 18, 40]. 
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Remark 7. At the moment we cannot present convergence rates for the Lavrentiev 
regularization under variational source conditions with the variant of a discrepancy 
principle, for which the convergence was shown in Section 2. Perhaps this is due to 
an intrinsic deficit of discrepancy principles in Lavrentiev regularization. Thus, we 
have considered the balancing principle (33) as an alternative a posteriori rule for 
proving rates, see Theorem 5. 


4. Convergence rates under approximate sonrce conditions. We assume 
that F locally satisfies an approximate range invariance condition, which means that 
there exist a bounded monotone (accretive) linear operator A G £(X,X) satisfying 
the condition 

(Ax,x}>0 for all xGX (37) 

and a ball Bp(x^) C 'D(F) of positive, but maybe small, radius p > 0 and a uniform 
constant 0 < c < 1, such that for all points x G Bp(x^) there are bounded linear 
operators M(x),N{x) G £(X,X) satisfying the conditions 

F(x) — F(x^) = AM(x)(x — x^) + X(x) and 
||M(:r)-J||<c<l, ||A^(x)||<a(||x-xt||) 


with (7 to be specified below (cf. (46)). If X does not vanish identically, we addi¬ 
tionally assume that 

Vx, X G Bp(x^), ||M(a;) — M{x)\\ < Lm ||a; — a;|| and AM{x) accretive 
- S|| < (1 - c)p, c-|-pLm<1 

holds. Note that by the convergence result from Section 2, there exists ^ > 0 such 
that for all 5 G (0,5] the regularized solution a;a.(<5) Bp(x^). 

Accretivity (37) of A has the consequence that (see, e.g.. Proposition 2.1 in [37]) 

||(A -I-a/)“^|| < — and ||(A-|- a/)“^ A|| < 1 for all a > 0. (40) 

a 

If the monotone operator F is differentiable at x G with Frechet derivative 

F'(x), then it is well-known that A := F'(x) is accretive and thus satisfies conditions 
(40). 


4.1. General convergence rates. The benchmark source condition for an accre¬ 
tive linear operator A G £{X,X) we consider here is 

x^ — X = Aw, for some w G X. (41) 

Moreover, we can use for arbitrary i? > 0 an approximate source condition of the 
form 

x^ - X = Awn + rn wr, xr G X, |]r(;fl|] = R, ||rfl|] = d(R), (42) 

with a distance function 

(i(i?) = min{ |]a;'^ — a; — Aui|| : ||u>|l < i?}. (43) 

If the benchmark source condition (41) fails and hence ||ri{|| > 0 for all i? > 0 in 

(42) , then the method of approximate source conditions using the distance function 

(43) applies. This method was developed for linear ill-posed operator equations in 
[19] (see also [11, 21] and ideas in [7]). An extension to nonlinear operator equation 
can be found in [17] and moreover in [10]. It yields convergence rates whenever 
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d{R), i? > 0, is continuous, positive, convex and strictly decreasing to zero as 
i? —>■ oo. This is the case if 

TZ{Af = X. (44) 

(see, e.g., [10, Lemma 3.2]). For accretive A this is equivalent to the injectivity of 
A due to TZ{A) 0 A/’(A) = X (cf. [16, Prop. 2.1.1 h]). 

Under the regularity conditions (41), (42), we will use the a priori parameter 
choice a ~ ^“^(<5) in the sense that there exists a sufficiently large constant r > 1 
such that 


6<T(j){a), and d(x ^(a)) < rc?(x ^(0 ^(<5)) if (41) is violated 
i(5 < < tS if (41) is satisfied 


(45) 


and the assumption 


a{e) < Co-ee ^(e) 

holds for some constant > 0 in (38), with the auxiliary functions 


(46) 


e{a) 


</>(a) = 

xiR) 


ad{x~^{a)) 
d{R) 


+ T)d{x ^(a)) if (41) is violated 
+ T)a if (41) is satisfied. 


(47) 

(48) 

(49) 


Note that d is monotonically decreasing by definition and x is strictly decreasing; 
if the benchmark source condition (41) is violated, then d is strictly decreasing and 
hence e is invertible and a is strictly monotonically increasing. Therewith, we arrive 
at the following rates result. 


Theorem 8. Let F : X ^ X be a hemicontinuous and monotone operator satisfying 
(38) with (37), (44) and (46), and, if N ^0, additionally (39). Then with the choice 
(45) we have 


llxl - xfj 


f 0{d{x ^(<5)))) */ (41) is violated, 

I 0{'/6) if (41) is satisfied. 

In particular, if F is Gateaux differentiable with locally Lipschitz continuous deriv¬ 
ative 

||F'(a;) — F'(a:’i)|| < T ||x — xfll for all xGBp{x^) (50) 

and (41) holds with A = F'{x^), then ||x^ ~ xiH = 0{VS). 


Proof. We first consider the case TV = 0 in (38), where the estimates are somewhat 
simpler. We rewrite (4) as 

xi-x^ 

= - (A + al)~^ {y - - F{x^) - A{xl^ - xf) + a{x^ - x)) 

= - (71 + al)-\y - /) - (71 + aI)-^A{Mixi) - I){xi - x^) 

— a{A + aI)~^Awji + a{A + al)~^r]i , 

with wr = w and = 0 in case of (41). Taking into account the assumption (38) 
(with N = 0) and the inequalities (40), we obtain the estimate 

Ika “ ^^11 ^ -\-aR-\- d{R), 

a 
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i.e., 

+ + (51) 

Considering the case that (41) is violated and equilibrating the terms containing R 
we obtain the best R in dependence of a: 

R = X-\a), (52) 

where x is defined as in (48). Equilibrating now the two terms in the right hand 
side of the resulting estimate 

\\xi - a;i|| < + 2d(x”^(a))) 

1 — c \a / 

yields the best a priori choice (45) (where <j) is defined as in (47)) and thus, the rate 
estimate 

\\xi-x'<\\=0{dix~\(l)~^{S)))) as (5^0. (53) 

If (41) holds, the estimate simplifies to 

< Y^(^ + “ll^ll) ^ ||w||)v^ 

where we have used (45). 

If N does not vanish identically, we (as in [37] ) introduce an intermediate quantity 
z solving the equation 

[AM{z) + al){z — x^) + a{x^ — a;) = 0 (54) 

with M as in (38). Solvability of (54) follows from Banach’s Fixed Point Theo¬ 
rem, applied to the reformulation 

z = 4>(z) = x^ + {A + al)~^ (^(-f ~ M{z)){z — x^) -I- a{x — a;^)^ (55) 

of (54). Indeed, 4> is a self-mapping on Bp{x^), since for any z G Bp{x^) 

||<I>(z) — a:^l| = ||(A -I- aI)~^A{I — M{z)){z — x^) -I- a{A + al)~^{x — x^)|| 

< ||(/-M(z)(z-xt)||+ ||x-xt|| 

<cp+ ||x-x'''||p< p, 

due to (40) and (39). Contractivity of $ follows from the estimate 

||$(z) - $(5)|| = \\{A + aI)-^A({I - M{z)){z - x^) -{I- M{~z)){'z - x^)) || 

< 11(7 - M{z)){z - 5)11 + ||(M(z) - M(5))(5 - xt)|| 

< (c-7Tm/o)||-2-5|| 


for any z, 5 G Sp(xl), where c -|- LmP < 1 by (39). 

Using (54), we can rewrite (4) as 

F(xf,) - F{z) + y -y^ + F{z) - F{x^) - AM{z){z - x'^') -I- a{xf^ 

By testing it with x„ — z, by monotonicity of F and (38), we obtain 


-^11 < - 
a 


a(\\z - x'l'l 


z) = 0. 


a 


(56) 
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Here, by employing (40) and the approximate source condition (42) with 
\\M{z)~^wr\\ < ||rfl;|| < d[R), we can estimate as follows, 


\\z — x^W = a\\{AM{z) + al) ^(a;’^—ai)|| 

R 

= a\\{AM{z) + aI)~^{AM{z)M{z)~^WR + r_R)|| < -- a + d[R), 

1 — c 


Using this and the triangle inequality in (56) we get 


Iki -x^ll < - + 


-(X 4 “ dyRj 4 - — 


a 1 — c a 

since a is monotonically increasing. By choosing the interdependence of i?, a, and 
5 in the best way as above, cf. (52), (45), we can achieve that 

d R 

a 4- d{R) < e{a) 


1 — c 


with e according to (49). So in order to maintain the rate from the case N = 0, 
we assume the function a : R"*" —>■ M"'" to be monotonically increasing and satisfying 
(46). This gives the estimate 

||xi - at II < e(a) + < (i + C^)eia), 

a 

which yields the claimed rates due to (45). □ 


Remark 9. If (41) fails, it follows due to (44) that the function d{R), R> 0, given 
by (43) is continuous, positive, convex and strictly decreasing to zero as i? —>■ oo, 
the rate function d{x~^{4'~^{S)) is strictly increasing for d > 0 and satisfies the 
limit conditions lim c?(x“^((()“^(i5))) = 0 and lim = 0. The latter 

condition shows that we have a lower convergence rate if the exact source condition 
(57) fails. 

This approach also works if d{R) stands for a majorant function to the distance 
function (43) which is strictly decreasing to zero as i? —>■ oo. 

If d is differentiable, then instead of equilibrating terms one can as well minimize 
the right hand side of (51) with respect to a and R, which leads to the same result 
at least in the Holder case considered in Section 4.3. 


4.2. Comparison to previous results. If the benchmark source condition 

— X = Aw, w G X, Ijwll = Ro > 0, (57) 

is valid, then (42) holds with WRg = w and d{R) = 0, Rq < R < oo. In that case 
we have the benchmark rate 

Ika “ as (5-)> 0, (58) 

whenever a ^ v^. This was already shown in [37] for the Frechet derivative F'{x^) G 
jC(X, X) as linear operator A, under the weaker nonlinearity condition (50) that we 
also use in Theorem 8. Note that for accretive A G C{X, X), condition (57) requires 
that — a; is orthogonal to the nullspace M{A). 

The range invariance condition (38) with N = 0 (cf., e.g., [24, 26, 34, 39]) is 
really a strong nonlinearity condition, but it has the advantage that A need not be 
exactly equal to F'{x^), but can be chosen from a wider variety of linear operators. 
For example, the nonlinearity condition involving constants ko,p> 0 such that for 
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all x,x £ Sp(x’t') C 'D{F) and h £ X there exist elements k{x,x,h) £ X with the 
property 

[F'(x) — F'(x)] h = F'(x) k(x,x,h) and ||fc(a;, x, < fco||x — a;|| (59) 

is frequently used (see, e.g., [34, conditions (10), (11)])- Due to the mean value 
theorem, this implies the range invariance condition (38) with A = F'(x^) provided 
that the smallness condition kg ||F'(a;^)|| < 1 holds. If the following inequality is 
satisfied for some constants 0 < k < 1 and fci > 0 , 

\\k{x,x^,h)\\ < ki ||a;-a;'''||'^||h||, 

then the smallness condition is fulfilled whenever p is sufficiently small. 

Condition (59) with fixed x := x^ occurs as Assumption A3 in [37]. Note that 
the operators A^, 0 < p < 1, mentioned at (5) are fractional powers of the accretive 
(sectorial) operator A, defined by a Dunford integral as 

OO 

APy ■- sin(p7r) y (60) 

0 

This approach was used there to prove Holder convergence rates 

llxf, — = O as 5 —>■ 0 , (61) 

with A := F'(x^) for Lavrentiev regularization and with an a priori parameter 
choice a ~ under the source condition 

x"^ — X = A^ w, w £ X, (62) 

(see also [33, Chapt. 8 ] or [31, Chapt. 1]). We will obtain this rate result in Sec¬ 
tion 4.3 under the range invariance condition (38) with a completely different proof. 

If (59) holds for all center elements x £ Bp(x^) and all x £ Bp(x'^), then one even 
has a range invariance condition (38) with A = F'(x) for all such x whenever 

\\k{x,x,h)\\ < fci jjx — fjl'' ||/i||, Wh £ X 

and p is sufficiently small. 

Extension of this type of conditions to (38) with N 0 allows for situations in 
which the ranges of the linearizations do not coincide exactly. Exact coincidence of 
the ranges can be quite restrictive if these ranges are non closed (as relevant in ill- 
posed problems). In the situation of dense but non-closed range, a small correction 
N might often be sufficient to bridge the gap between ranges of linearizations at 
different points. 

4.3. More on Holder convergence rates. 

Proposition 10. Let F : X ^ X be a hemicontinuous and monotone operator 
and suppose that a range invariance condition (38) holds for some accretive linear 
operator A satisfying (44) (with additionally (39) if N ^ OJ. Provided that the 
benchmark source condition (57) fails for all Rq > 0, we have a Holder convergence 
rate (61) for Lavrentiev regularized solutions for an a priori parameter choice a ~ 
( 5 p+t ij there is some source element w £ X such that the fractional power source 
condition (62) with A^ from (60) is satisfied with 0 < p < 1. 
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Proof. Under the assumptions of the proposition, estimate (53) is applicable and 
thus we only have to show that, for all 0 < p < 1 and under (62), 

^(x~^('/'~^(^))) = 0((5^) as (5 —0. (63) 

In this context, we exploit the formulas 

\\{AA* + XI)-^AP\\<^ and \\A*{AA* + XI)-^Ap\\ < (64) 

2 A— 

which are valid for all A > 0 and constants Ci,C 2 > 0 depending on p S (0,1). 
These formulas are a consequence of the moment inequality for sectorial operators 
A (cf. [16, Prop. 6.6.4]). 

The minimum problem (43) for verifying d{R) can be reformulated by the La¬ 
grange multiplier method, where for any A > 0 the element v\ € X is the uniquely 
determined minimizer to 

||x^ — X — Av\\'^ -I- A||i;|p —?► min, subject to v £ X, 
which can be written as 

vx = (AIM + XI)-^A*{x^ -x)= A*{AA* + XI)-\x^ - x). 

Since by (44) the range of A is dense in X and since the source condition (57) 
fails, we have that the strictly decreasing function 0(A) := ||ua||^ = ||A*(Ayl* -|- 
A/)“^(a;’^ — a;) 11 ^ satishes the limit conditions 

lim 0(A) = -|-oo and lim 0(A) = 0. 

A—^0 A—^“|”Oo 

Then there is a uniquely determined Ar > 0 such that 0 (Ar) = and d{R) = 
||x'f — X — AuarII. Based on (62) and on the first inequality in (64), we derive the 
estimate 

d{R) < \\Xr{AA* + XrI)-^APw\\ < ArIKaIaI* + Ar/)-Mp||||u;|| < CiA|||u;||. 
We claim that the following decay rate can be established 

diR)<KR^, 0<R<oo, (65) 

for some constant K > 0. Then by using the majorant function K Rp^ as d{R) in 
formula (53), one obtains the required condition (63) with x(i?) ^ i?p^, 4>{a) ^ 
and under the a priori parameter choice a ~ (5p+t. Thus, it remains to show 
(65). By the second inequality in (64) we have 0(A) < UlllwlpA^”^ and thus the 

value Xmaj = (cTpJll) ^ ^ > 0 solving the equation C' 2 ||u>||A^^j = R satisfies the 
inequality Ar < Xmaj- This implies the estimate (65) and completes the proof. □ 

Remark 11. Due to the identity 

TZ(AP) = TZ((AA*}P''^) 

(cf., e.g., [32, Lemma 1]), which is valid for all 0 < p < 1, the source condition (62) 
is equivalent to 

x'^ — X = (AA*)P‘^^ w, w £ X, (66) 

and Proposition 10 remains true if (62) is replaced with ( 66 ). In that version the 
proof can be found as a direct consequence of Theorem 3.2 in [11] (extendable also 
to non-compact operators A) by using A instead of A*, while taking into account 
that A £ C(X,X). This also yields the decay rate (65) and hence the required 
Holder rate result. 
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Remark 12. In order to provide a relation between the variational source con¬ 
ditions from Subsection 3.2 as well as from Subsection 3.1 and Holder source 
conditions of range-type, we consider the special case of a monotone operator 
F := A G C{X, X) and we set for simplicity x := 0. 

Firstly, let A = A* he a, self-adjoint operator and let the range-type source 
condition (62) hold for some 0 < p < 1 , which implies the convergence rate (61). 
Note that only for 0 < p < 1/2 we have the inequality chain 

= {w,A^x) < ||r(;|| ||^^a;|| 

based on the interpolation inequality. By using Young’s inequality 
ab < ^ + ^ with a = \\x\\'^~'^p, b = llmll (Ax,x)p, ^ and p = this 

yields (for all 0 < p < 5 ) the variational source condition 

{x\x) < + Q + 11^11 ^ {Ax,x)^ (67) 

of type (28) occurring in Subsection 3.2 with p = and valid for all x G X. 

Evidently, from (67) we derive in the sense of Subsection 3.3 the Holder convergence 
rate 

||x^(^) — a;^|| = 0((5^) = 0(i5^) if q;((5) ~ S ' =^p+t^ ( 68 ) 

which coincides with (61) for all 0 < p < 1 / 2 . 

Secondly, consider an accretive A which is not necessarily self-adjoint, and let 
for some 0 < p < 1 the source condition 

x^ = {A*)P w, w G X, (69) 

hold. According to [16, Proposition 7.0.1 (e)], one has {A*)p = {A^)* and the 
inequality chain 

= {{A*Yw,x) = {w,APx) 

for some c > 0 , which is a consequence of the moment resp. interpolation inequality 
for monotone operators, see [31, Corollary 1.1.19] or [16, Proposition 6.6.4]. Then 
Young’s inequality yields, for all 0 < p < 1, the variational source condition 

(a;^x) < ll^^f + (■^) Ikll^ \\Ax\\^ 

of type (18) occurring in Subsection 3.1 with ^/(t) = and valid for all x G X. 
This gives for 0 < p < 1 

l| 2 ^a( 5 ) — a;^ II = 0(5 ) if a{6) ^ . ( 70 ) 

From [32, Theorem 1 and Lemma 1] it can be seen that under the source condi¬ 
tion (69) the obtained rate (70) is not optimal for 0 < p < 1/2 and can be improved 
to as in ( 68 ), because Ti{AP) = TZ{{A*)p) for 0 < p < 5 . On the other 

hand, for 1/2 < p < 1 the ranges TZ{Ap) and TZ{{A*)p) are in general different for 
non-self-adjoint accretive operators A as the Volterra operator shows (see Subsec¬ 
tion 5.1). From [32] one can also see that for 1/2 < p < 1 and under the source 
condition (69) the rate (70) is only optimal in the special case p = 1. 


5. Examples. 
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5.1. Applying Lavrentiev-type variational source condition to linear prob¬ 
lems with Volterra operator. This subsection provides an example where a 
Lavrentiev-type variational source condition does provide the desired convergence 
rate. We consider A = L^(0,1) and the linear Volterra operator (simple integration 
operator) A : X —>■ X, defined by 

[Aa;](s) := f x(t)dt, 0 < s < 1 , 

Jo 

which is accretive (monotone), but not self-adjoint. For our study let us focus on 
the specific solution = 1. This violates the source condition = A*v required 
for the rate ||x^ — x'^|| = O ^(5^^ in [28], as for this purpose x^ would have to vanish 

at the right boundary point. Moreover, x^ = 1 does not fulfill x^ = A^v which 
would imply the same estimate along the lines of [37]. In fact, we have that x^ = 1 
satisfies x^ G TZ{A^) = TZ{{A*Y) for any p G (0,1/2), but we have x^ ^ TZ{A^) 
and x^ ^ 'TZ[[A*Y) for p G [1/2,1]. This can be derived from the explicit structure 
of the ranges of the fractional powers of the Volterra operator given in [15], and 
one can simply see that the corresponding ranges TZ{Ap) and TZ{{A*Y) are different 
for all 1/2 < p < 1. Indeed, according to [16, p.234. Prop. 8.5.5] we can rewrite 
= A^v as the Abel-type integral equation 

x^(s) = (A^u)(s) = } . [ (s — tY~^vit) dt, 0 < s < 1 . 

r(p) Jo 

Due to the identity r(p)r(l — p) = the well-known explicit solution formula 

for the Abel integral equation attains the form 


which yields v(t) = ^ for x'^ = 1 so that u G L^(0,1) iff p < 

However, a variational source condition (28) resp. (29) is satisfied, as shown 
below. 

Let z(s) := [Ax](s), 0 < s < 1. Integration by parts yields 
z(t)x(t)dt = f z{t)dz{t) 

Jo 


(Ax,x) = / 
Jo 


This implies for all x G A 


{z(l)f- [\{t)dz{t) = {z{l)Y-{Ax,x). 
Jo 


and hence we have for x^ = 1 and all x G A 


(Ax,x) = i^y x(t)dt^ 

1 and all X G A 
(x^,x) = f x{t) dt < V2{Ax,x)^, 

Jo 


thus (28) holds with /3i =0, P 2 = and p = ^. Then Theorem 5 implies the 
convergence rate 

Ika ~ =0 a si. 

The same result is also a consequence of Theorem 8 using approximate source condi¬ 
tions. Namely, we have d{R) ~ for the distance function 
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d{R) = inf{||x'l' — Aw|| : ||it;|| < R}, which immediately follows from the argu¬ 
ments of [13, Example 5] by replacing there the function •\/2cos(i — l/2)7rt) with 
•\/2sin(i —l/2)7rt) in order to execute the transfer from A* to A in the eigenfunctions 
of the singular system. Then we have from Theorem 8 

\\xi -x^\=0 = O ((5^) . 

5.2. A parameter identification problem in an elliptic PDE satisfying the 
local range invariance condition. Consider identification of the source term q 
in the elliptic boundary value problem 

— Au + £{u) = q in n 

(71) 

u = 0 on dn 


from measurements of u in fl, where ^ : R —^ R is some Lipschitz continuously 
differentiable monotonically increasing function and fl C R^ a smooth domain. 
The corresponding forward operator F : L?{yi) —>■ C g i—>• u, is indeed 

monotone, since 


{F{q)-F{q),q-q) 


/ {u — u){q — q) dx 

Ja 


j (u — u) ^—A(u — 




dx 


= + (C(m) -C(m))(m- w)) 

> ||V(m- > 0 


dx 


(72) 


(which by the way still does not imply strong monotonicity of E, since this would 
require a lower bound in terms of H? ~ 911^2 (q))- Lipschitz continuity of F follows 
from the fact that w = F{q) — E(gl) can be written as the solution of 


—Aw -f ^qW = q — q^ in 
w = 0 on dn 


(73) 


where 


iq= [ +t(E(g) - u'f))(it > 0 

Jo 


(74) 


with = F{q^). Indeed, by testing with w, integration by parts and using Poincare- 
Friedrichs’ as well as Cauchy-Schwarz inequalities we get 

Il'“'lli2(n) A ||Vw|||2(n) 


< 


= / w{q-q^)dx <111x11 L^n)\\q-q^lmn) ■ 
Jn 


dx 


According to Theorem 1.4.6 in [2], the operator F is actually maximal monotone, 
since it is monotone, continuous (hence hemicontinuous) and its domain is the space 
L‘^{Q). We define the linear operator A : L^(n) —>■ L‘^{Q) as h i—>• Ah = v with v 
solving 

—Au + £'(u^)v = h in n 
n 

z; = 0 on ai Z 


where solves (71) with q = q^. 
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We claim that F satisfies the range invariance condition (38) with N = 0. In 
order to show this, we use the representation oi w = F{q) — (73), (74) from 

above, as well as (75), which yields that F{q) — F{q''') = AM{q){q — gi) holds with 

M(g) : ^ , M{q) = (-A + ^'(u^) id)(-A + ^ id)"! . 

Note that for any a € a > 0 a.e., the linear operator 

(—A + a id) : 77^(17) fi Hq{D,) —>■ L^{Vl) , u i—>■ —Am + au 

is well-defined and continuously invertible and both A and AM{q) are accretive 
since for any a G i^(n), a > 0 a.e., (so e.g., a = or a = |g), h G and 

V := (—A -I- a id)~^h, we have 


((—A-f a id) ^h,h) = f {\\7v\'^ + av^) dx > 0. 

Jn 

The difference of M{q) to the identity can be estimated as follows: 

\\Miq) - III = II ((-A + r(«^) id) - (-A + 4 id)) (-A + 4 id)-i I 
= ll(C'(«^)-|,) id(-A + 4id)-i|| 

ll(r(«n-4)(-A + 4id)-VlU= 


= sup 

0//GL2(n) 


\L^ 


sup 

o^wem{n)nH^in) || - Am; -I- |,m;||l2 
< Ga sup - — - 


where we have used elliptic regularity, i.e., the above mentioned mapping properties 
of (—A -I- a id), and continuity of the embedding i7^(fl) —>• L°°{Q) with norm 
■ Here one has 


lir(u^)-1,11^2 = 


< 


n 


< 


^'(m^) — ^'(m^-|- t(M — M^)) dt) dx 
L|/ y (y t{u — u^)dt^ dx 

y (y dt J {u — u^Ydt^dx 


L^, 


3 “ 3 

for some constant Lj/ > 0. Consequently, we obtain 


"1?-^ A£ki|9-«'lG, 


\u - u’ I|l2 < 


||M(g) - /|| < CACg.^^. ||g - gt |1^. < c 


for all g G Hp(g'i), where c = is smaller than one, provided 

that p is sufficiently small. 
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We now verify the Lavrentiev specific source condition (29) for the special case 
of Holder rates ip{t) ^ G [0, ^]. More precisely, we will show that 

(which in case /r < |; is equivalent to q^ — q G hence does not impose 

any boundary conditions on q^ — q) is sufficient for this variational source condi¬ 
tion, provided the domain is sufficiently smooth to admit the elliptic regularity and 
embedding bounds 

=: < oo , ||(-A)^“^|| j/ 2 - 2 M(Q)_>,i 2 (Q) =: (7^2^ < > 

ll(~A) =: < oo , ||id||ii-- 2 ;.(n)_,,i 2 (n) =: < oo , 

where — A is the Laplace operator equipped with homogeneous Dirichlet boundary 
conditions. For this purpose we use the fact that by (72), for any q G Ai := Bp{q"^) 
we have 

7 := {F{q)- F{q'^),q- q'^) > || V(m - u 1 ')||^ 2 (n) > Ik - 

as well as 


ii^(u) - ?(k)ii^-2.(o) < c% ii^(u) - au^)h^n) = j^iau) - au^)r dx 

< - ^{u^)){u - wt) dx < 

The first of these two estimates can be further made use of in the interpolation 
estimate 

II A / + Ml ((~A)(m'I' — u),v) 

II - A(u^ - u)||j^^ 2 m(q) = sup 


= sup 

^GC-(n)\{o} 


^)GC^(a)\{o} lkllff=»‘(n) 
((—A)^“''(u'l' — u), {—A)'^v) 


< C'^C'^2/i Ik^ - M|k2-2;.(n) 

^ C'^C'^2/i Ik^ ~ “lliL2(n) Ik^ ~ '“lliiqn) 

< C^,Ct2,{C^ II - A(ut - ix)|U2(n))^-^k(l + Cpf)^Y , 

where we can further estimate 


II - A(Y - u)llL2(a) = Ilk -q - (kk) - kk)lk2(n) < P + \kkV 7 - 
This altogether yields 

(k - 9, k - 9 ) = ( 9 ^ - 9, -A(k -u) + - ^H) 

^ Ik^ - 9llH„"'^(n) (ll - - u)lk-=»^(n) + Ikk) - kk)|k-2;.(f2)) 

< Ik'^' — 

■ {cilCt2^C!^(P + V^k7))'-"k(l + Crf)7)'' + 

=: ifY) = 0 ( 7 ''). 
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6. Conclusions and further work. In this paper, we have shown convergence 
rates for Lavrentiev’s regularization method under variational and approximate 
source conditions for linear and nonlinear inverse problems with monotone forward 
operators. In particular, we have proposed a new variational source condition that 
seems to be quite appropriate for the Lavrentiev setting. 

To compare the capability of the different source conditions, we present in the 
following table the best possible rates in the linear case. 



condition 

rate 

(a) 

= Aw for some w G X 

0{5^) 

(b) 

x^ = A*w for some w G X 

0{5^) 

(c) 

{x^,x) < /?! |a; 1^ -1- I32{Ax,x)^ for all x G X 

0((53) 


As the example from Section 5.1 shows, (c) implies neither (a) nor (b). However, 
(c) is implied by a fractional source condition = A^w in case of a self-adjoint 
operator A, cf. Remark 12. The question of rates beyond those stated above 
appears to be a challenging one, which we intend to investigate further on. 

Note that no additional restriction on the nonlinearity of F is needed as regards 
variational source conditions- not even differentiability - as already observed in [22] 
for Tikhonov regularization. On the other hand, we had to impose some local range 
invariance condition in order to prove rates with approximate source conditions. 

Lavrentiev’s method in Banach spaces will be subject of further research in light 
of the few aspects considered in this Hilbertian setting. 
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